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We consider N vicious walkers moving in one dimension in a one-body potential v(x). Using the 
backward Fokker-Planck equation we derive exact results for the asymptotic form of the survival 
probability Q(x, t) of vicious walkers initially located at (2:1, . . . , xjsr) — x, when v(x) is an arbitrary 
attractive potential. Explicit results are given for a square- well potential with absorbing or reflecting 
boundary conditions at the walls, and for a harmonic potential with an absorbing or reflecting 
boundary at the origin and the walkers starting on the positive half line. By mapping the problem 
of N vicious walkers in zero potential onto the harmonic potential problem, we rederive the results 
by Fisher [J. Stat. Phys. 34, 667 (1984)] and Krattenthaler et al. [J. Phys. A 33, 8835 (2000)] 
respectively for vicious walkers on an infinite line and on a semi-infinite line with an absorbing wall 
at the origin. This mapping also gives a new result for vicious walkers on a semi-infinite line with 

jV(jV-l) 

a reflecting boundary at the origin: Q(x, £) ~ t 2 



I. INTRODUCTION 

Consider N symmetric random walkers which annihi- 
late on meeting each other but otherwise do not interact. 
This concept of short-ranged interacting random walkers 
was introduced by M. E. Fisher as the vicious walkers 
model Q. 

One of the main properties of interest in this model 
is the survival probability Q(x, t) that none of the 
N vicious walkers with initial position coordinates 
(xi,X2, . . • ,xn) — x has met another up to time t, i.e. 
none of them has been annihilated up to time t. 

Fisher and Huse 0, determined the survival prob- 
ability for N vicious walkers moving on an infinite line. 
For large times, Q(x, t) decays as a power: 



Q(x,t) ~ r 



(1) 



Interesting results also arise when further conditions are 
imposed on the movement of the vicious walkers by the 
use of absorbing or reflecting walls, where all walkers 
are initially located on the same side of the boundary 
(the case where there are walkers on both sides decou- 
ples into two independent problems). While Fisher [| 
found the survival probability problem of vicious walk- 
ers with an absorbing boundary at the origin only for 
N = 2, Krattenthaler et al. Q were able to determine 
exact asymptotic forms for N vicious walkers starting 
from equi-spaced lattice points: 



Q(x,t) ~r 



(2) 



By evaluating the scaling limit Katori and Tanemura |j] 
showed that this asymptotic behaviour holds for arbi- 
trary initial conditions on a continuous line. 

In this paper we introduce the interesting problem of N 
vicious walkers moving in an attractive one-body poten- 
tial v(x), i.e. the full potential function has the separable 
form V(x) = Y^ii=i v i x i)- Treating both time and space 
as continuous, we investigate the survival probability of 



N vicious walkers with equal diffusion constants D. The 
equation of motion for walker i is taken to be 



dV 

dxi 



(3) 



where the Langevin noise r\{ (t) is a Gaussian white noise 
with mean zero and correlator 



(Vi{t)Vj(t'))=2DS ij 6(t-t'). 



(4) 



For a square- well potential of width L we consider three 
different combinations of absorbing and reflecting walls 
and find an exponential decay for the survival probability 
of the general form Q(x, £) ~ e~ 0Nt . For two reflecting 
walls the exponent 9n is determined to be 



qRR 



= D 



7T 2 N(N -1)(2N -1) 
1? 6 ' 



(5) 



In the case of one reflecting and one absorbing wall we 
obtain 



N 



D 



N(2N +1)(2N -1) 



L 2 



12 



(6) 



while for two absorbing walls the exponent of the asymp- 
totic decay is: 



q aa 



= D 



7T 2 N(N + 1)(2N + 1) 
L2 6 ' 



(7) 



An interesting potential which turns out to a powerful 
tool is the problem of TV vicious walkers in the harmonic 
potential V(x) = fx 2 . The asymptotic behaviour for 
large times is determined to be an exponential decay in- 
dependent of the diffusion constant: 



JV(JV-l) 
* n = a. 



(8) 



This result also provides a mechanism to determine the 
survival probability of N vicious walkers on an infinite 
line in a simple way. By mapping the zero-potential 
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problem to the harmonic potential problem, we derive 
Fisher's result [lj and also the result by Krattenthalcr 
et al Q with an absorbing wall at the origin. Further- 
more we are able to obtain, to our knowledge, a new 
result for the survival probability of N vicious walkers 
on a semi-infinite line with a reflecting boundary at the 
origin, which decays as: 



Q(x,t) ~f 



N(N-l) 



(9) 



The outline of the paper is as follows. In section II we 
present the method for a general one-body potential v(x), 
while in section III we give explicit results for square- well 
and harmonic potentials. In section IV we revisit the case 
of zero potential, obtaining the known results, and a new 
result for a system with a reflecting boundary, through a 
transformation to the harmonic problem. Section V is a 
short conclusion. 



II. THE METHOD 

The dynamics of a random walker, with position co- 
ordinate Xi, moving in a potential V(x) is described, in 
continuous space and time, by the Langevin equation @ 
with noise correlator Q. 

The probability <5(x, t) that all N vicious walkers, 
i = l,...,N, have survived up to time t, given that 
they started at {xi}, satisfies the corresponding backward 
Fokker-Planck equation: 



dQ(x,t) 

at 



N 

i— 1 1 



y, dT/(x) 3Q(x,t) 
f-f dxi dxi 



(10) 

For convenience we start by defining the survival proba- 
bility q(xi, t) of just one random walker moving in a po- 
tential restricted by the imposed boundary conditions. 
This survival probability q(xi,t) satisfies the backward 
Fokker-Planck equation 



dq(xj,t) 
dt 



d 2 



dv(xi) dq(xi, t) 



dxi 



dxi 



(11) 



where we have used the relation V(x) = Y^i v ( x i) for a 
one-body potential. For any such potential, the back- 
ward Fokker-Planck equation is separable in time 
and space. Let us call these separable solutions, i.e. the 
solutions of Eq. satisfying the relevant boundary 

conditions, single-walker basis functions. They have the 
form qj(xi,t) — Uj{xi) exp(— Xjt), where A,, is the decay 
rate associated with basis function j, and these rates are 
ordered such that Ai < A2 < A3 . . .. 

For N non-interacting walkers moving in the same 
potential, the JV-walker basis functions for the survival 
probability take the form of products of N single-walker 
functions, each with a different space variable x^. Since, 
however, we are investigating vicious walkers the mu- 
tual annihilation property must be respected. Since two 



walkers die when arriving at the same cc-coordinate, the 
boundary condition Q(x\, .. .,x n ,t) — when Xi = Xj 
for any i ^ j must be respected. This property is en- 
sured by constructing Q(x, t) using antisymmetric com- 
binations of products of N single- walker functions, anal- 
ogous to the antisymmetric construction of the wavcfunc- 
tion of fermions [l| . The N- walker basis functions of the 
vicious walker problem with N walkers have, therefore, 
the form 

Q^'-'^x,*) = .lei .1' ' N . (12) 
where the elements of the N x N matrix A are given by 



(13) 



The full solution Q(x, t) is a linear superposition of these 
basis functions with coefficients determined by the initial 
condition. 

To solve the problem of N vicious walkers in an ar- 
bitrary potential, therefore, we need to find the single- 
walker basis functions qj (&, , t) appropriate to the im- 
posed boundary conditions. For an absorbing boundary 
at x = a the functions qj (xi , i) must satisfy 



q 3 (xi = a,t) = 0. 



(14) 



For a reflecting boundary at x = b the boundary condi- 
tion for the backward Fokker-Planck equation is in gen- 
eral 



VQ(x,t)-n = 0, 



(15) 



where n is normal to the reflecting boundary 0. In the 
one-dimensional case, this expression implies 



dqj_ 

dxi 



0. 



Xi—b 



(16) 



Clearly these boundary conditions are also satisfied by 
the functions Q n ''"' w (x, t), since the latter is just an 
antisymmetrised product of single- walker basis functions. 

Consider now the late-time limit, t — > 00. Each 
antisymmetrised product in the expression for Q(x, t) 
contains N different relaxation factors exp(— Xjt). The 
slowest-decaying term in the sum, therefore, is the term 
in which the relaxation rates are Ai, A2, . . . , A at. It fol- 
lows that, asymptotically, 



Q(x, t) oc detB 1 - 2 '-'^ exp(-6 N t), 



(17) 



where B is just the N x N matrix with elements 

Bum — u n {x m ) (n,m — 1, . . . ,N), i.e. it is constructed 
using the N slowest-decaying single-walker basis func- 
tions, and the total decay rate is 



7jV 



N 



(18) 



The following sections provide some applications of this 
general result. 
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III. RESULTS FOR VICIOUS WALKERS IN A 
POTENTIAL 

In this section we discuss two examples of N vicious 
walkers in a potential and determine the decay of the 
survival probability Q(x, t). 



A. The square-well potential 

Consider a square-well potential which has two walls 
of infinite potential, one at the origin and the other at 
x = L, and vanishes between the walls. A vicious walker 
restricted to move between the walls satisfies the back- 
ward Fokker-Planck equation: 



dq(xi,t) _ D d 2 q{x l} t) 



dt 



dx 2 



(19) 



2. One reflecting and one absorbing wall 

For an absorbing wall at the origin and a reflecting 
wall at x = L the boundary conditions are satisfied by a 
Fourier sine series with basis functions 



q„{xi,t) = exp - 



(2n + l) 2 ir 2 Dt 
4L 2 

x sin (^-(271+1)^) , n = 0, 1, . . .(24) 

Analogous to the preceding case the survival probability 
for all TV vicious walkers is constructed and the asymp- 
totic survival probability for large time is evaluated us- 
ing E J =o 1 ( 2i + !) 2 = N ( 2N + !)( 2Ar - !)/ 3 t0 S ive the 
asymptotic decay 



Q(x, t) <~ exp 



ir 2 Dt N(2N+ 1)(2N- 1) 
~ 2 12 



(25) 



This equation can be solved in general by separation of 
variables, which amounts in this case to writing the solu- 
tion as a spatial Fourier series. Different solutions result 
from the various sets of boundary conditions imposed by 
the property of the walls. 



3. Two absorbing walls 

In the case of two absorbing walls the basis functions 
have to vanish at both at x = and at x — L. A Fourier 
sine series is therefore appropriate, with basis functions 



1. Two reflecting walls 

For two reflecting walls the spatial derivative of q(xi, t) 
must be zero at x = and x = L. In this case, there- 
fore q(xi,t) is given by Fourier cosine series with basis 
functions 



, n 2 ir 2 Dt 
q n {x t ,t) = exp [ — — | cos 



n = 0, 1, 



(20) 

The survival probability is constructed as a superposition 
of antisymmetrised products of these basis functions: 



Q(x,t) = j2---J2 cil "' iN detAil "" 



l N 



£•••££" ex p 

ii in 

X del ir '\ 



ir 2 Dt 
L 2 



N 

71=1 



where 



DJl,...,ljV 



COS 



(-■ \ 



(21) 



(22) 



To evaluate the long-time behaviour we keep only the TV 
longest-lived modes, given by the N smallest values, i = 

0, 1, . . . , N - 1 of *„. Using Etq 1 ? = N ( N - !)( 2iV " 
l)/6 we obtain, for the asymptotic time-dependence, 



Q(x,t) — exp ( - 



ii 2 Dt N(N — 1)(2N- 1) 
~L 2 6 



(23) 



n 2 n 2 Dt 



sin (J^-nXij , n = 1, 2, . 



q n (xi,t) = exp 

\ ±j- i \jj / 

(26) 

This is very similar to the result for two reflecting bound- 
aries, except that the spatial functions are sines so the 
sum begins with n = l. The large-time behaviour Q(x, t) 
is given by 



Q(x.,t) ~ exp ( - 



ir 2 Dt N(N+ 1)(27V + 1) 
~L 2 6^ 



(27) 



Before proceeding to the harmonic potential, we note 
that the inequalities 2iV(JV - 1)(2JV - 1) < N(2N + 
1)(2N - 1) < 2iV(iV + 1)(2N + 1), for all N > 1, im- 
ply that for a well of given size the decay is fastest with 
two absorbing boundaries and slowest with two reflecting 
boundaries, as is intuitively clear. 



B. The harmonic potential 

A harmonic potential V(x) — , . 
which the backward Fokker-Planck equation for the one- 
walker basis function reads 



dq(xut) d 2 



a Xl ^^. (28) 
oxi 



This equation can be transformed into an imaginary- 
time Schrodinger equation by the substitution q(xi,t) = 
exp(ax 2 /AD)ij)(xi,t) to give 



AD 



il){x u t). (29) 
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This equation has solutions of the form tp(xi,t) = 
e~ u(xi), where u(xi) satisfies the ordinary differential 
equation 



D- 



dxi 



AD 



u(xi) = -Xu(xi). 



(30) 



This equation is equivalent to the time-independent 
Schrodinger equation for the harmonic oscillator. The 
eigenvalues and eigenfunctions of this eigenvalue prob- 
lem are well known: see for example a similar problem in 
ref. 6]. The eigenfunctions have the form 




a \ a 2 \ 

2D W-4D^)> 



(31) 



where the functions H n (x) are the Hcrmitc polynomials 
defined by 



Hn(y) = (-iTe* 



(32) 



The corresponding eigenvalues are A n = na, where n = 
0,1,2,.... The original basis functions q(xi, t) are given 
by q n (xi,t) = H n (a^Vsj) exp(-X n t). 

Applying the antisymmetrisation process to determine 
the survival probability of N vicious walkers in a har- 
monic potential we obtain the asymptotic time depen- 
dence: 



N—l 



Q(x, t) ~ exp I — at i 



giving 



<9(x, i) ~ exp ( —at 



N(N- 1) 



(33) 



(34) 



This approach is readily extended to the case where 
there is a reflecting or absorbing boundary at x — and 
all the walkers start on the same side of the boundary (if 
there are walkers on both sides, the problem decouples 
into two independent problems). For a reflecting bound- 
ary, the boundary condition u'(0) = selects only the 
even-numbered Hermite polynomials, n = 0, 2, 4, . . ., and 

/ N-l \ 

Q(x, i) ~ exp I — at 2i J 

= exp[-atN(N -1)} (reflecting wall). (35) 

For an absorbing boundary, the boundary condition 
u(0) = selects the odd-numbered Hermite polynomi- 
als to give 

Q(x,i) - exp ^-ai^(2i - l)j 

= exp [~atN 2 ] (absorbing wall). (36) 

In the following section we show how these results can 
be used to compute the survival probability of N vicious 
walkers in zero potential, with and without an absorbing 
or reflecting wall, by mapping the problem back to the 
oscillator problem. 



IV. VICIOUS WALKERS ON A LINE 

Here the case of N vicious walkers restricted by no 
potential is investigated. This problem can be solved in 
a quite simple way by mapping it to the problem of N 
vicious walkers in a harmonic potential and using the 
previous results. Again, we consider the Langevin equa- 
tion ©, but with V(x) = 0, and let all N vicious walkers 
start to move at time t = to. We introduce the following 
mapping from x, t to the new coordinates X,T by 0: 



t = t Q e 1 



Then the Langevin equation J3Jl transforms to 



d -^§ 1 = -\x i {T) + UT) 



(37) 



(38) 



where &(T) = V^A^ e T / 2 m {t e T ) is a Gaussian white 
noise with mean zero and correlator 



(Zi(T)Z j (T'))=5 ij 6(T-T'). 



(39) 



The corresponding backward Fokker-Planck equation in 
the new coordinates is 



1 N 

IE 



2 ^ dX 

i=l 



2 ^ 



dXi 



dQ(X,T) 
dT 

i — i - i—i 

(40) 

where the space coordinates are now the starting points 
of the vicious walkers, given by: 



Xi(T = 0) 



Xijtp) 

V2DTo 



(41) 



In the new coordinates this problem looks identical to 
the harmonic potential problem with a = 1/2 and D = 
1/2. Hence the asymptotic (in time) solution for the 
survival probability of N vicious walkers is, according to 
our previous results, 



T 



N-l 



Q(X,T) ~ exp i-^Yl A^ B " 



(42) 



i=0 



where (B H ) nm = H n -i(X m fy/2) and n,m = 1,...,N. 
Mapping back to the original coordinates (x, t) leads to 
the survival probability 



Q(x,t) 



2 i=0 ' 



det^ 



(43) 



where (B L ) nm = H n _ 1 (x m /2 y / Dt ), with n,m = 
1,..,,N. The determinant deti? L is proportional 
to the Vandermonde determinant : det B L — 



(Dt ) 



-JV(JV-l)/2 



n 



°3\> 



and all ^-dependence 



drops out, as it must, to give the long-time behaviour 



Q(x,t) ~i" 



(44) 
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which is just the result Fisher obtained Qj. But our ap- 
proach to the problem also gives a simple way to ob- 
tain expressions for the survival probability for N vicious 
walkers with an absorbing or reflecting wall at the origin 
(and all walkers starting on one side of the wall). 

The essential arguments have been given in the preced- 
ing subsection. For an absorbing (reflecting) boundary, 
only the odd (even) basis functions contribute. Note first 
that the Fisher result l|44|) follows immediately from i|34|) 
on setting a = 1/2 and T = ln(t/to). The detailed discus- 
sion above was given mainly to show how the arbitrary 
scale to drops out. To obtain the asymptotic results for 
a reflecting or absorbing wall at the origin, we can sim- 
ple make the same replacements in Eqs. PB*|) and ifHri)) 
respectively. For the absorbing boundary, we recover the 
result of Krattenthaler et al. [3j : 

Ml 

Q(x,t)~t 2 (absorbing wall), (45) 
while for a reflecting wall we obtain 

N(N-l) 

Q(x,i)~t 5 — (reflecting wall). (46) 

The latter is, to our knowledge, a new result. 

As a final comment we note that the case where the 
absorbing or reflecting wall moves, with a displacement 
x w = ct 1 / 2 , is also amenable in principle to exact anal- 
ysis. The change of variable (|37f) maps the problem to 
one where the N walkers move in a harmonic oscilla- 
tor potential, and the absorbing or reflecting wall is at 



a fixed position in the new coordinates. This problem 
has been analysed for a single walker and the sur- 
vival probability decays as t~ e , where the exponent 9 
is found to vary continuously with the amplitude, c, of 
the wall displacement. The same qualitative features 
will be present for N vicious walkers. For a reflecting 
(R) or absorbing (A) boundary, one will obtain a de- 
cay exponent Or. a = N(N — l)/4 + (c, N), where 
fR,A(—°°,N) = 0, corresponding to a rapidly receding 
wall, which will be equivalent to no wall at all, and 
f R (0,N) = N(N - l)/4, /a(0, N) = N(N + l)/4 cor- 
respond to a static wall. 



V. CONCLUSION 

In this paper we have derived the exact asymptotics 
for the survival probability of vicious walkers moving in 
a square well potential and a harmonic potential with 
various combinations of absorbing and reflecting walls. 
The results for a harmonic potential have been used to 
find the properties of free vicious walkers (zero potential) 
through a change of variables, and a new result obtained 
for the case of a single reflecting boundary. Compar- 
ing all results for each potential one recognises that the 
survival probability decays faster when the number of 
walls is increased, with absorbing walls causing a faster 
decrease than reflecting walls, in accord with intuitive 
expectations. 
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